High fidelity all-optical control of quantum dot spins: 
detailed study of the adiabatic approach 
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Confined electron spins are preferred candidates for embodying quantum information in the solid 
state. A popular idea is the use of optical excitation to achieve the "best of both worlds", i.e. 
marrying the long spin decoherence times with rapid gating. Here we study an all-optical adiabatic 
approach to generating single qubit phase gates. We find that such a gate can be extremely robust 
against the combined effect of all principal sources of decoherence, with an achievable fidelity of 
0.999 even at finite temperature. Crucially this performance can be obtained with only a small 
time cost: the adiabatic gate duration is within about an order of magnitude of a simple dynamic 
implementation. An experimental verification of these predictions is immediately feasible with only 
modest resources. 



I. INTRODUCTION 

Electron spins in quantum dots (QDs) are promising 
candidates for quantum computation (QC) due to their 
long intrinsic decoherence tiniest Optical control of such 
spins, via auxiliary exciton (electron-hole) states, offers 
the promise of fast gating times. The 'figure of merit', i.e. 
the ratio of decoherence time to gate time, would then 
be extremely high. This exciting possibility has been 
discussed in several proposals recently^ 3 -^^ 

However, direct exploitation of the excitonic degree of 
freedom, which we term dynamic optical control, may ad- 
versely affect the spin coherence: During the gate opera- 
tion the quantum information is partially carried by the 
excitons, which are subject to aggressive decoherence. It 
has been suggested that adiabatic control could avoid this 
problem, by ensuring that the qubit remains encoded in 
low-lying states throughout the processi^ In this paper 
we present calculations for both forms of gate under the 
combined effect of the principal decoherence mechanisms, 
photon emission from excitonic recombination and acous- 
tic phonon interaction. In contrast to previous work^ 
we derive a full master equation (ME) solution for the 
gate dynamics that simultaneously incorporates all deco- 
herence channels. 

Our principle interest is in evaluating the performance 
of the adiabatic approach, and the Markovian master 
equation technique that we develop is tailored for this 
purpose. We will show that, although the adiabatic- 
ity condition introduces the possibility of Landau-Zcncr 
transitions, this docs not significantly impair gate perfor- 
mance. Error rates less than 10 -3 can be achieved even in 
the presence of multiple decoherence channels. In order 
to provide a meaningful context for these results, we will 
also consider the performance of a simple dynamic ap- 
proach in the same parameter regime i.e. relatively weak 
driving. However this form of gate is not the primary 
focus of our study, and we will neglect more advanced 
dynamic models and techniques such as pulse shaping, 8 
and also the effects of pure dephasing that may become 
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FIG. 1: (Color online) Pauli blocking effect: heavy holes oc- 
cupy the lowest energy valence band states with spin m z = 
±3/2 whereas light holes have spin m z = ±1/2 at a slightly 
higher energy. Therefore, a + light cannot excite electron hole 
pairs if the qubit is in the state |0) with m z — —1/2. On the 
other hand, if the qubit is in the 1 1) configuration the excita- 
tion is possible leading to a three particle trion \X) state. 



significant in the case of ultrafast addressing.— 

Our predictions for both the dynamic and the adiabatic 
gates can be tested immediately; experiments would only 
involve a single QD and one laserJ^ Such a demon- 
stration would provide an important step towards im- 
plementing more general gates based on adiabatic ap- 
proaches 



II. MODEL 

Consider a self-assembled QD that is doped such that 
one excess electron, the spin of which is the qubit, per- 
manently occupies the lowest energy state of the conduc- 
tion band. If the dot is irradiated with a + polarized laser 
light only one of the electron spin configurations is com- 
patible with the creation of an additional exciton due to 
the Pauli blocking effect^ (illustrated in FigUJ. The re- 
sulting three particle state is called a trion and denoted 
\X). In this case, the logical |0) state (defined as spin 
down | j)) is unaffected by the laser pulse. Exploiting 
the selective coupling |1)(= | f)) to \X) thus allows for 
exciton-mediated spin manipulation. 
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In a frame rotating with the laser frequency oji, and 
after making a rotating wave approximation (RWA), 
we obtain the Hamiltonian of such a QD in the basis 
{|0),|1),|X)}(S=1) 

F s = A|X)(X| + |(|l}(X|+h.c), (1) 

where h.c. denotes the Hcrmitian conjugate, A = luq — 
lji is the detuning between uji and the exciton creation 
energy u>q, and f2 is the laser- QD coupling strength. 

Hamiltonian {TJ) provides a valid description of the 
driven QD, regardless of whether a classical laser field 
or a quantum mechanical laser mode in a coherent state 
is assumed^ However, the eigenstates of the joint sys- 
tem in a fully quantum mechanical treatment contain the 
photon number N of the laser mode 2 ^, giving rise to the 
following representation in the so-called 'dressed basis'^ 

\-) N = cos0|l,iV+ 1) -nm9\X,N), (2) 
\+) N = Bm0\l,N + l) + con0\X,N), (3) 

6 = \ arctan (§) ■ ( 4 ) 

The states |— )jv and |+)jv are conveniently grouped and 
referred to as a manifold M.(N). In this picture, there is 
a ladder of manifolds, each manifold separated from its 
neighbours by the energy of a laser photon uji J£ Hence- 
forth, we shall label the |— }jv state simply as |— ) when- 
ever information about the photon number in the laser 
mode is not needed. 

Based on Hamiltonian ([TJ) , a Rz{4>) gate (|0) — ► |0), 
|1) — ► e^|l)) can be performed in two different ways. 
First, a resonant Rabi flop allows one to achieve any angle 
of rotation provided control of the driving laser phase 
is possible^ Second, slow switching of an off-resonance 
laser beam can be used to achieve adiabatic following of 
instantaneous system eigenstates. As depicted in Fig [2] 
population initally in |1) follows the |— ) state as 9 goes 
from zero to some value 9 max and returns back to |1) 
as 9 goes back to zero. During this process, the energy 
difference between |— ) and |0) , 

E m = i (A - VA 2 + tt 2 ) , (5) 

causes an accummulation of phase of |1) relative to |0). 
This kind of adiabatic following can be achieved by using 
a laser pulse with slowly changing Gaussian field ampli- 
tude f2(i) = flo cxp[— (i/r) 2 ] with a constant detuning 
A. 



III. RADIATIVE DECAY 

Excitonic lifetimes up to a nanosecond have been re- 
ported.— Therefore, the problem of spontaneous photon 
emission has previously often been assumed to be in- 
significant compared with other decoherence channelsi^ 
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FIG. 2: (Color online) Energy of the dressed states: 

shown as a function of the detuning A in arbitrary units. 
The |— ) state tends to |1) and |+) to \X) for a large positive 
detuning. At resonance, A = 0, the dressed states are an 
equal superposition of the bare basis states. The decoupled 
|0) state is also shown, its energy remaining constant. Adia- 
batic following: population in the |0) and |1) states follows 
the respective instantaneous eigenstates if the mixing angle 
changes sufficiently slowly. |0) is unaffected by the laser pulse 
but |1) goes into |— } and back. Due to the energy shift of the 
|— ) state, a dynamical phase accumulates relative to |0). 



However, to achieve adiabaticity, our quantum gate must 
be performed slowly and the finite excitonic lifetime be- 
comes a relevant factor in limiting performance. Spon- 
taneous emission then causes qubit dephasing, although 
it does not cause qubit relaxation since |0) is always un- 
coupled to the photon emission process. 

Dynamic gates also suffer from such decoherence since 
any population in the excited state is susceptible to ra- 
diative decay at the rate of the inverse natural lifetime 
Tq. Let S denote the overall 'number of expected de- 
cays' for a square 27r pulse that is defined by the integral 
of the excitonic population over the pulse duration. For 
undamped, resonant Rabi oscillations, the population of 
the excited state obeys (|A)(A|) = sin 2 (f2/2i), so that 

S = r o y sin 2 (fii/2)di = -r . (6) 

In principle a stronger driving reduces S. However, in- 
creasing the pulse amplitude is technically demanding 
and can lead to population leakage from the trion sub- 
space. 

In the adiabatic scheme, radiative transition rates be- 
tween dressed states [Eqs. © and ([3])] are required to 
describe spontaneous photon emission processes. These 
are obtained analogously to the case of an atom inter- 
acting with a laser pulsed As can be seen in Fig [31 the 
total decay rate T m from |— ) into the adjacent manifold 
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FIG. 3: Radiative decay from manifold A4(N) to M(N — 1). Left: the uncoupled basis. The energy between the two states in 
each manifold is the detuning A and the solid arrows correspond to absorption and stimulated emission processes, whereas the 
wavy arrows denote spontaneous emission. Right: Allowed spontaneous emission transitions between the dressed states. The 
energetic splitting in each manifold is the effective Rabi frequency fi' = V A 2 + fi 2 and the spacing between adjacent manifolds 
is the laser frequency uji. The left emission process is at frequency lui + Q', the two centre-lines emit at the frequency of the 
laser and the right at coi — Q'; this is the Mollow triplet well known in quantum optics.— 



is given by the sum of two processes |— )jv — * \— )n-i arL d 

\-) N - \+) N -i^ 

- 2 4 cos 2 6 + T Q sin 4 9 = T Q sin 2 6. (7) 



r m = r n sin 



Obviously, Y m decreases as 9 gets smaller. On the other 
hand, a smaller 9 entails a prolonged gating time. By 
using Eq. ^ and performing a series expansion in fi/A 
around f2/A = 0, we obtain 
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A similar expansion for the energy shift E„ 
state yields 



of the 
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(9) 



The time required for a Rz(tt) operation is r = Tr/E m . 
The expected number of decays during the gate is simply 
the product of time and decay rate, yielding 



3 = s r 



o- 



(10) 



Hence H decreases with increasing A at the cost a longer 
gating duration r. We verify this result in Fig. [H 
which compares the result Eq. (jTUJ) with a full numeri- 
cal simulation, showing an excellent agreement whenever 
A > 217. We have checked our results with an appro- 
priate Markovian ME equation 2 ^, which reconfirms that 
further detuning the laser reduces radiative decay in the 
adiabatic scheme, while it takes a stronger coupling CI to 
achieve the same for the dynamic gate. 



IV. PHONON INTERACTIONS 

Interactions with vibrational modes of the surrounding 
lattice cause twofold decoherence. First, they dephasc 
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FIG. 4: (Color online) Any significant population of the ex- 
citonic state makes decay events inevitable. The expected 
number of decays (for To = 0.01 ps _1 ) during an adiabatic 
Rz(it) gate is shown as a function of A. We depict the ana- 
lytical prediction (red) along with full numerical solutions, for 
which fio is given in meV. The inset shows the corresponding 
pulse durations r in ps. 



the qubit in a similar way to radiative decay and second 
they can lead to gate failure through state relaxation. 
We focus on deformation potential coupling to acoustic 
phonons as the main phonon induced dephasing mecha- 
nism for slow excitonic processes ; 16 ^ 7 since it dominates 
over the much weaker piezoelectric coupling ] i 18 

In addition to the electron spin qubit, the trion state 
\X) also consists of a strongly bound exciton. Therefore, 
its interaction with phonons is given b y 16 ' 17 



# ep = |X)(X|]T 5q (a q + h.c.): 



(11) 



where a q is the annihilation operator for a phonon 
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with wavevector q and g q is an effective excitonic cou- 
pling strength ,g q = (M°V[ip e (r)] - M%P[i) h {r)])M. The 

Mq h = D e / h |q| y/h/ '(2/iV"w q ) are the deformation po- 
tential coupling strengths for electrons and holes^ with 
/i being the mass density, V the lattice volume, and 
D e / h the respective electron and hole coupling constants. 
V[4> e / h ] denotes the form factor of the electron / hole 
wavefunction. 

We proceed by transforming Eq. pip into the di- 
agonal basis of Hs (Eq. ([T])) and by then writing it 
in the interaction picture with respect to both Hs and 



q ^q" q Uq 



yielding, 



us' ,q 



— iuj't 



h.c 



3q 



h.c.) ,(12) 



where J G {0, A}, P = cos 2 6>|+) (+| + sin 2 6\-)(-\ and 
P A = — sin0cos0|— )(+|. A = V^ 2 + A 2 is the energy 
difference between the dressed states. 

The ME is derived in the usual wayi£ by integrating 
the von Neumann equation for the density matrix g of the 
joint system and tracing over the phonon modes. This 
results in an integro-differcntial equation for the qubit 
density matrix p: 



(13) 



The Born-Markov approximation is now performed, 
which relics on two assumptions. First, there is no back- 
action from the small system on the much larger bath, 
meaning the joint density matrix factorizes at all times 
Q = p® pb- Second, the bath relaxation is assumed to be 
rapid and so we replace p(t') by p(t), thereby neglecting 
any memory effects^ Further, we assume that system 
dynamics occurs on a timcscale much faster than the de- 
coherence processes, J(A) <C A, allowing us to perform a 
RWA and leading to a ME in Lindblad forim 19 i 23 

p = J(A) ([N(A) + l]D[P A ]p + N(A)D[Pl] P ) . (14) 

D[L]p = Lptf - l/2(tf Lp + pL^L) is the 'dissipator' of 
the ME, N(u>) = (exp(cj/fc B T) - 1) _1 describes the ther- 
mal occupation of the phonon modes, and J(u>) is the 
phonon spectral density: 



J(u) = 2ir} y |gq| r5(w-w q ). 



(15) 



The Lindblad operator Pq has been dropped because 
the spectral density vanishes for u> = 0. Consequently, 
the ME in the Born-Markov approximation provides a 
suitable description of phonon-assisted transitions but 
does not describe pure dcphasin g l 7 '"^ in the diagonal ba- 
sis. Pure dephasing is an intrinisically non-Markovian 
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FIG. 5: (Color online) Phonons affect dynamic and adiabatic 
quantum operations differently. To characterize decoherence, 
we present the purity of the system's density matrix during 
a R z (it) operation for deformation potential coupling at T = 
5 K. The dynamic gate uses Q = 0.1 meV and time is shown 
in units of 2n/Q,. Adiabatic gates are performed with Qo = 
1 meV and the units of A are meV. We have scaled the time 
axis by dividing by 6r and adding 0.5 (thus mapping the 
actual integration interval of — 3r to 3r to the scale of this 
figure). The inset shows the deformation potential spectral 
density of GaAs [Eq. ([TBJ]. 



process^ that can be interpreted as a partial 'which-way- 
measurement' of the qubit by the environment, thus caus- 
ing decoherence^ However, it can be adiabatically elim- 
inated if system dynamics proceed no faster than at a 
characteristic timescale of the order of 1 ps^^i 9 - 

The calculation of the spectral density J(oj) re- 
quires a microscopic model for the carrier wave- 
functions. The precise wavefunction shape is not 
important 2 ^ so we proceed by using ground state solu- 
tions to a harmonic confinement potential of strength 
162 meV at 2.5 run from the centre of the dotiSI 
^Vfc(r) = (4//!^)"^ exp(-r 2 /(2d 2 //i )), for which 

d e /h = {ft/ yf™^Jhc) 1/2 and c = 8.3 x 10~ 3 J/m 2 . For 
GaAs>ii±£ D e = 14.6 cV, D h = 4.8 eV, p = 5.3 g/cm 3 
and c s = 4.8 x 10 5 cm/s. The effective electron and hole 
masses are m e = 0.067 mo, rrih — 0.34 mo. Assuming 



a linear phonon dispersion w q 



|q|, the spectral den- 



sity then takes a super-Ohmic form (shown in the inset 
of Fig. © 



J(w) 



2pc 5 s 



D 



-FT* * )• ( 16 ) 



The exponential cut-off terms on the RHS of Eq. JTJJ]) 
are related to the finite size of the QD and filter out 
high frequency phonons with wavelengths too short to 
interact with the dot. The cut-off frequencies are: uj e ^h = 
VZcs/de^h and uj eh = 2c s / ^/dfTdf . 
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By examining the structure of Eqs. (fT4f and (fl6|) we 
conclude that the influence of phonon-induced decoher- 
ence becomes small as A approaches zero and can be ex- 
ponentially suppressed for values of A beyond the spec- 
tral cut-off Lu e / h . Outside the ultrafast regime only the 
first route of decreasing A via ft lends itself to the dy- 
namic scheme but this is contrary to the requirement of 
avoiding radiative decay. On the other hand, operation 
beyond the cut-off is realizable for an adiabatic gate with 
a large enough detuning. Using the purity of the system's 
density matrix, tr[p 2 ], to characterize system decoher- 
ence, Fig. [5] shows that phonon decoherence can indeed 
be overcome with the adiabatic approach and that in gen- 
eral performance improves as the detuning increases 2 ^. 



V. LANDAU-ZENER TRANSITIONS 

Landau-Zener (LZ) transitions are non-adiabatic tran- 
sitions between eigenstates approaching an anticross- 
ingj22 they are a source of error in the adiabatic scheme 
since they cause population to transfer to excited states, 
subjecting it to fast decoherence during and after the 
gate. LZ transitions can be suppressed by using longer 
gating times. Contrary to the adiabaticity conditions de- 
rived in Refs. [H and [j, which are only valid for a linear 
sweep through resonance, we take a more general ap- 
proach. 

To derive an adiabaticity condition consider the 2LS 
Hamiltonian ([I]). The transformation to the basis of in- 
stantaneous eigenstates, 



U{6) 



cos V — sin ( 
sin 9 cos 



(17) 



is time-dependent. The Hamiltonian transforms as H = 
WHU + i (£f7 + ) U, yielding 

H = A"|-)(-| + A+l+X+l + + h.c.) , (18) 



where A ± = 1/2(A ± y/A 2 + ft 2 ) denote the instanta- 
neous eigenenergies of the states |— ) and |+). The off- 
diagonal terms give rise to LZ transitions, which are 
therefore suppressed if 



\0\ < IA+-A- 



(19) 



Using Eq. (01 for 9, we obtain a general condition for 
adiabaticity: 



ftA -ftA 
2(A 2 + ft 2 ) 3 /2 



< 1. 



(20) 



For n = r2o e x p[~(V T ) 2 ] an d constant A, we deduce that 
inequality (|2"D|) can be satisfied by demanding 



ft /A 2 < r. 



(21) 



Adiabatic following is therefore always achieved in the 
limit f! < A, where r oc A over a broad range of param- 
eter space. 
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FIG. 6: (Color online) Comparison of the dynamic versus the 
adiabatic gate for a Rz(t) operation by showing the overall 
gate fidelity. Left: fidelity of the dynamic operation as a 
function of the coupling strength Q. Right: fidelity of the 
adiabatic operation for a fixed f2o = 1 meV as a function of 
the detuning. 



VI. OVERALL GATE FIDELITY 

In order to bring together the results so far, we now 
calculate the fidelity of the Rz{n) gate with a combined 
ME describing both spontaneous photon emission and 
phonon-induced processes. The gate is performed on the 
input state = (|0) + |l))/\/2 and ideally produces 
|V>-> = (|0) - |1))/V2 as its output state. The fidelity is 
defined as T = (ip-\p\ip-), where p is the density matrix 
of the system after the gate has finished. Again, we use 
the material parameters of GaAs and assume a radiative 
decay rate of r = 0.01 ps _1 . 

The left panel of Fig. [6] shows the fidelity of the dy- 
namic gate as a function of ft for different temperatures. 
Towards small values of ft, the fidelity is limited by the fi- 
nite excitonic lifetime, whereas phonon-induced processes 
become important for larger values of ft. The highest fi- 
delity here is below 0.95 at absolute zero and decreases 
at finite temperatures. For much higher values of ft be- 
yond the phonon spectral cutoff (~ 10 meV), we would 
expect the fidelity to improve rapidly, neglecting the ef- 
fects of pure dephasing. However, as outlined above, our 
approach is not suited to studying the regime of ultra- 
fast driving, and we refrain from showing actual values 
at large ft. 

The fidelity of the adiabatic gate as a function of A 
for a fixed value of fto = 1 meV is shown in the right 
panel of Fig. [5] Towards larger values of A, T increases 
and asymptotically approaches unity as radiative decay 
is more effectively suppressed. Phonon-induced processes 
arc temperature dependent and predominantly occur for 
small values of A. The intermediate peak visible at 
higher temperatures is related to the dip in the spectral 
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density of the deformation potential (see inset of Fig. somewhat lower gate fidelity. 



VII. CONCLUSIONS 

We have performed a comprehensive Markovian deco- 
herence study of an exciton-mediated spin phase gate in 
quantum dots. Within an adiabatic scheme all princi- 
pal sources of decoherence can be effectively suppressed 
in the same limit of weak off-resonant laser driving, and 
a phase gate fidelity of 0.999 or better may be possi- 
ble even at finite temperature. In contrast, for similar 
driving amplitudes, a dynamic gate with typical operat- 
ing parameters suffers strong decoherence that leads to 
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